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Abstract 

We present an explicit polynomial formula for evaluating the principal logarithm of 
all matrices lying on the line segment {1(1 — t) + At : t G [0,1]} joining the identity 
matrix / (at t = 0) to any real matrix A (at t = 1) having no eigenvalues on the closed 
negative real axis. This extends to the matrix logarithm the well known Putzer's method 
for evaluating the matrix exponential. 
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1 Introduction 

Given a nonsingular matrix A G IR nxn , any solution of the matrix equation e x = A, where e x 
denotes the exponential of the matrix X, is called logarithm of A. In general, a nonsingular real 
matrix may have an infinite number of real and complex logarithms. If A has no eigenvalues 
on the closed negative real axis then A has a unique real logarithm with eigenvalues in the 
open strip {z G C : — it < Imz < tt} of the complex plane (see, for instance, [5]). This unique 
logarithm may be written as a polynomial in A and is called the principal logarithm of A. It 
will be denoted by log A. 

The problem of computing the principal matrix logarithm has received some attention 
in recent years (see, for instance, pQ, [2], 0, |~D and 0). In part, this interest has been 
motivated by the applications of the matrix logarithm in areas such as Systems Theory and 
Control Theory. The above cited papers list some applications. 
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As far as we know, most of the methods proposed for computing the principal logarithm 
are approximation methods . Unlike the matrix exponential case, for which several closed 
forms based on polynomial representations have been studied (see, for instance, [Zj, jH], 
and [10J), little attention has been paid to closed forms for the matrix logarithm. 

In this paper, we find for the matrix logarithm the analogue of the well known Putzer's 
method jU] for evaluating the matrix exponential. Assuming that for t G IR the spectrum of 
I — At does not intersect IRq , we consider the curve t i— ► log(J — At) in IR raxn . Using the 
coefficients of a polynomial p(X) of degree k such that p(A) = 0, every matrix in that curve 
will be written as a linear combination of the matrices I,A,---, A k ~ 1 ) in the following way: 

log(/ - At) = h{t)I + f 2 (t)A + ■■■ + f k {t)A k -\ 

where the coefficients fx, - ■ ■ , fk are integrals of certain rational functions. 

We find this simple method suitable for teaching purposes because the topics required for 
understanding it (basically, eigenvalues of matrices and integration of rational functions) are 
usually taught in the first years of undergraduate courses. We recall that, in contrast, other 
methods proposed for evaluating the matrix logarithm require advanced theory, such as Schur 
decompositions, matrix square roots and matrix Pade approximants. 



2 A polynomial formula for the matrix logarithm 

Given A G IR nxn , let p(X) = A fc + ciA fc " 1 + - • - + Ck-iX + Ck be a polynomial with real coefficients 
such that p(A) = and let 

" • • • -c fc 

Ik-l '■ 
-ci 

where I m denotes the m x m identity matrix, be the companion matrix of p(\). Examples of 
polynomials p(X) such that p(A) = are the characteristic polynomial of A (k — n) and the 
minimum polynomial of A (k < n). 

Before stating our main result, let us define the following subset of IR: 

V = [t G IR : a (I - At) n IR = 0} , 

where cr(X) stands for the spectrum of X and A is a given n x n matrix. For each t G IR, the 
eigenvalues of I — At are of the form 1 — Xt, with A G <j{A). Since non real eigenvalues of A 
always give rise to non real eigenvalues of I — At, it is enough to consider real eigenvalues of 
A to obtain a more clear description of the set T>. Thus, we may write 

V = {t g IR : 1 - At > 0, VA G a{A) n IR}. 
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Let Am = raax(a(A) (1 IR) and A m = mm(a(A) D IR). Assuming that A has both positive 
and negative real eigenvalues, we have V =]j-, ^-[. If A does not have negative eigenvalues 
then V =} — oo, ^-[ and if A does not have positive eigenvalues then V =}j-, +00 [. In any 
case, P is an open interval. 

Theorem 2.1 Suppose that the above notation holds and that the vector function [fi(t), ■ • ■ , fk(t)] T 
is the solution in T> of the initial value problem 

(I-Ct)x(t) = -e 2 , x(0) = 0, (1) 

where e 2 = [0 1 • • • 0] T . Then 

log (/ - At) = h{t)I + f 2 (t)A + ■■■ + f k (t)A k -\ (2) 

for all t G V. 

Proof. The function X(t) = log (I - At) is differentiate for all i G V and X(t) = -A(I - 
At)^ 1 (see (6.6.14) and (6.6.19) in [3]). Besides, X(t) is the unique solution in V of the initial 
value problem 

(I-At)Y(t) = -A, Y(0) = 0, (3) 

where Y{t) G IR nxn . 

Let [hit), • • • , f k (t)f be the solution of flU) in V and define P(i) := f\{t)I + f 2 {t)A+ ■ ■ ■ + 
fk(t)A k ^ 1 . In the following we show that P(t) is also a solution of (JBJ). Clearly P(0) = since 
/i(0) = 0, Vi = l,---,A;. 

Since the vector function [/1 • • • /fc] T satisfies (JTJ), a little calculation lead us to the system 

A + c k tf k = 

~tfl + f 2 + Ck-ltfk = -1 

< -tf 2 + h + c k -2tfk = . (4) 

. -t/fc-i + (1 + cit)f k = 
Using the equations of and the identity A k = —c\A k ~ x — • ■ • — Ck~\A — c^J, which follows 
from the Cayley-Hamilton theorem, we may write 

(I-At)P(t) = (I-At)(f 1 I + f 2 A + --- + f k A k ~ 1 ) 

= fJ + (h ~ tf x )A + ■-. + (/,- tf k -i)A k ^ - tf k A k 
= (fi + c k tf k )I+(-tf 1 + f 2 + c k -itf k )A + --- 

■■■ + (-tfk- 1 + (l + c 1 t)f k )A k - 1 
= -A 

Since (JHJ) has a unique solution, it follows that P(t) = log (I — At), m 

Since the coefficients functions in (J2J) are solutions of (fl|). we can obtain formulae for 
fi, i = 1, • • • , k, by solving the first equation for fi and substituting it into the second 
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equation, solving the second equation for f 2 and substituting it into the third equation and 
proceeding similarly until the last equation. The result is 



fi — —Ckfkt 

fi = ~t l 2 — fk Sj=l Ck-i+jft j i = 2, ' ' ' , k — 1, 
fl 



or, equivalently, 



k 



./': 



l+c 1 t+---+c k t k • 



1 1 + Cl t + - + C k t k 



f _ -t *-cif c k -it" — O U 

■> l ~ l+cit+-+c k t k ' ' ^' ■ " ' ^ 



/fc 



l+eit+---+c fc t fe ' 

We note that the constants arising in the integration process to find f\, % — 1, ■ ■ ■ , k, can 
be evaluated according to the identities /i(0) = 0, i = 1, ■ ■ ■ , k. 

We now summarize the previous discussion in the next corollary. 

Corollary 2.2 Given A e \R nxn , let p(X) = X k + CiA fc_1 H h c fe _iA + c k be a polynomial 

with real coefficients such that p(A) = 0. If V = {t G IR : a (I — At) fl IR = <f>}, then, for all 
t G V, 

log (/ - At) = h(t)I + f 2 {t) A + ■■■ + f k (t)A k -\ 
where f\, ■ ■ ■ , fk are differentiable functions in T> given by 

fl(t) = Jo l+ Cl C s +...+c fe s fe ds 

M) = si "'^:^ v "' ds ' i = 2,-,*-l (5) 

fk(t) = f 1+cls *... +CfeS fe ds. 

Remark 2.3 There exists a relationship between the polynomial p(X) and the polynomial 
g(A) = 1 + CiA + ■ • ■ + CfcA fc in the denominator of the functions under integral symbols in (jSJ): 
g(A)=A fc p(l/A). 

Remark 2.4 The indefinite integrals in (JHJ) may be obtained explicitly because we are dealing 
with rational functions. We note that many calculus textbooks provide methods for evaluating 
integrals of these kind of functions. Also, symbolic software packages like Mathematica, Maple 
or Derive are able to compute them. 
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Remark 2.5 For A such that cr(A) D IR = 0, equation (j2J) allows us to find an explicit 
formula for evaluating the logarithm of all matrices on the line segment joining I (at t — 0) 
to A (at t = 1): 

{I(l-t)+At:te [0,1]}. 

Indeed, 

log(I(l-t)+At) = log(I-(I-A)t) 

= f\(t)I + f 2 (t)(I — A) -\ h f k (t)(I - A) k -\ 

Obviously, this formula holds not only for all t G [0, 1], but also for any t such that a (I — (I — A)t)n 
IRo = 0- In particular, for t = 1 we may compute directly log A: 

log A = h(l)I + / 2 (1)(J — A) H h / fc (l)(I - Af-\ 

3 Example 

To illustrate the method proposed, we consider the matrix 





7 


4 


-4 


A = 


4 


7 


-4 




-1 


-1 
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To compute log A we have to work with the matrix I — A. The spectrum of / — A is 
{ — 11,-2,-2} and its minimum polynomial is p(A) = A 2 + 13A + 22. Applying directly 
(JSJ), we have 

/•' 22s ft —1 

= / 1 i iq i oo 2 rfs ' = / 1 i 1, i oo 2 

io 1 + 13s + 22s^ Jo 1 + 13s + 22s z 

Evaluating the integrals, we may write 

log (/-(/- A)t) = h(t)I + f 2 {t){I - A), 

where t £ [0, 1] and 

hit) = ^ln(l+2t)-|ln(l + m) 

Therefore 

log.4 = fi(l)I + Ml)(I - A) 

- (ln3 + ?ln(i))/ + Iln(i)(/-A). 

Acknowledgements: The author would like to thank Professor F. Silva Leite for sug- 
gesting this problem and for useful conversations on this topic. 
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